Modeling Equilibrium Clusters in Lysozyme Solutions 



Frederic Cardinaux*, Anna Stradner, Peter Schurtenberger 
Department of Physics, University of Fribourg, CH-1700 Fribourg, Switzerland. 



o 
o 

(N 

wo: 
< 



Francesco Sciortino, Emanuela Zaccarelli 
Dipartimento di Fisica and CNR-INFM-SOFT, Universita di Roma La Sapienza, 

(Dated: February 6, 2008) 



Italy. 



We present a combined experimental and numerical study of the equilibrium cluster formation in 
globular protein solutions under no-added salt conditions. We show that a cluster phase emerges 
as a result of a competition between a long-range screened Coulomb repulsion and a short-range 
attraction. A simple effective potential, in which only depth and width of the attractive part of the 
potential are optimized, accounts in a remarkable way for the wavevector dependence of the X-ray 
scattering structure factor. 
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Competition between short-range attraction and long 
range repulsion provides an efficient way to stabilize ag- 
gregates whose shape and characteristic size result from 
the delicate balance between these opposing forces S 
El [3 El Under appropriate external conditions, par- 
ticles interacting with such a mixed potential may form 
equilibrium cluster phases. This is a state of matter in 
which the stable structure of the solution is characterized 
by the presence of equilibrium aggregates of particles, the 
colloid analog of micelles 0. 

Cluster phases have been recently observed in colloidal 
systems as a result of the competition between short- 
range depletion attraction and long-ranged electrostatic 
repulsion^, 0]. Confocal microscopy has provided de- 
tailed information on the size and shape of these clusters. 
Theoretical and numerical studies suggest that cluster- 
phases can also be observed in different systems such as 
star-polymer solutions^] or charged lyposomes^. The 
typical signature of the equilibrium cluster phase is a 
pre-peak in the structure factor, signaling a preferential 
distance set by the competing forces on different length 
scales. Such a feature has been recently reported in solu- 
tions of globular proteins 0, 0], implying a generality of 
the mechanism by which bulk aggregation is disfavored 
and hnite-size clusters are formed and persist in equi- 
librium. This similarity suggests the possibility of an 
approach to globular proteins based on the assumption 
of an effective interaction potential, in which the compe- 
tition between repulsion and attraction is built in. 

In proteins, short-range attraction is normally at- 
tributed to a combination of van der Waals attraction, 

I — ~i I I 

hydration forces, and hydrophobic interactions yjj, II 1| . 
While in colloid-polymer mixtures such short-range at- 
traction is well characterized in terms of polymer size 
and concentration, in protein solutions it is poorly un- 
derstood. However, clear evidence for its presence is 
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provided by studies at increa sing ionic strength of the 
second virial coefficient 0, ll3L IT^j and by the de- 
termination of gas-liquid coexistence lines 0, 0, Il7| . 
which are metastable with respect to crystallization as 
for short-ranged attractive colloids fis| . Long-range re- 
pulsion arises from screened electrostatic interactions, as- 
sociated to the net charge of the protein in pH-controllcd 
protein solutions. 

Among globular proteins, lysozyme has become the 
prototype for scientific investigations. Under no-added 
salt condition, the lysozyme structure factor shows a 
clear cluster pre-peak whose position is essentially inde- 
pendent of protein concentration and weakly dependent 
on temperature^]. Hence, it is particularly important 
to assess under which conditions an effective potential 
can be designed which accounts for such typical features. 
Previously, theoretical studies have attempted to study 
lysozyme solutions for different solution parameters (e.g. 
pH, density, salt concentration), either with an effective 
continuum model for electrostatic interactions between 
proteins carryin g d iscrete charges [l^ or with explicit 
primitive models At high ionic strength, modeling 
for the calculation of the phase diagram usually relies 
on purely attractive potentials |ifj|. Some models have 
incorporated the 'patchyness' and/or the non-sphericity 
of the interactions [lfl EH ], associated with the discrete 
distribution of charged and hydrophobic sites, whose rel- 
ative distances could compete with the relevant distances 
incorporated in an effective spherical potential. But, to 
our knowledge, for none of these models, despite the large 
number of involved parameters, an equilibrium cluster 
phase has been predicted. 

In this Letter we present a combined experimental 
and numerical study on the equilibrium cluster forma- 
tion in suspensions of a globular protein in the absence 
of added salt. We propose an effective simple poten- 
tial for the lysozyme-lysozyme interaction in water based 
on a short-range attractive potential complemented by a 
Yukawa screened electrostatic repulsion. Parameters in 
the Yukawa potential are fixed by the known size and 
charge of the protein and the composition of the sol- 
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vent. In contrast to previous studies on colloid-polymer 
mixtures where the contribution of background ions is 
dominant [i^, HH , counterions coming from the surface 
charge induce a highly concentration-dependent behav- 
ior of the amplitude and screening length of the Yukawa 
potential, which is properly taken into account [24\- The 
depth and width of the potential is chosen by best-fit 
with the structure factor at one reference concentration, 
and kept fixed for all other studied state points. The re- 
sulting potential is capable of accurately describing the 
measured static properties at low and intermediate con- 
centrations. 

The globular protein used is hen egg white lysozyme 
(Fluka, L7651). Its molecular weight is 14.4kDa and 
its shape ellipsoidal with linear dimensions of 3 x 3 x 
4.5nm[l^. For simplicity, we neglect the asymmetry 
in shape, and we model it as a sphere of diameter 
(7 = 3.4nm. Lysozyme is dispersed in a solution of 
D 2 (99.9%, Cambridge Isotope Laboratories) contain- 
ing 20mM Hepes buffer salt. A detailed description of the 
sample preparation procedure can be found elsewhere 3 . 
The pH = 7.8 is adjusted with sodium hydroxide and 
held constant within ±0.1 units for all volume fractions 
investigated. Under these conditions, the net charge 
of the protein is known from titration experiments to 
be Zq = +8e 25]. The ionic strength of the solvent, 
estimated by conductimetry, is 8mM, corresponding to 
a Dcbye screening length at infinite protein dilution of 
£ ~ 3.4nm. SAXS experiments are performed with a 
pinhole camera (NanoSTAR from Bruker AXS) equipped 
with a sealed tube (Cu K Q ), a thermostatically regu- 
lated sample chamber and a two-dimensional gas detec- 
tor. Samples are measured at volume fractions ranging 
from cj) = 0.085 to 0.201, where (j) = 7rper 3 /6, with p the 
protein number density, is systematically measured by 
UV-visible spectroscopy. 

We perform molecular dynamics simulations (MD) of 
a system composed of N = 2500 particles of diameter 
a and mass m in a cubic box of size L, as a function 
of volume fraction <j> and temperature T. The excluded- 
volumc term plus the short-range attraction are modeled 
for simplicity with the generalized Lennard-Jones 2a — a 
potential 



V SR (r) = 4e 



(1) 



where a essentially controls the width of the attraction 
and the steepness of the hard-core, while — e is the depth 
of the potential. The parameters a and e are chosen 
as units of length and energy respectively. The short- 
range potential is complemented by a long-range repul- 
sion, modeled with a Yukawa potential. Its amplitude 
and screening length are fixed by the experimental condi- 
tions and follow the generalized one-component macroion 
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FIG. 1: Potential Vtotir) as a function of interparticle distance 
r for various cj>. The potential well depth is 4.22fc_eT and its 
range is « 3.6%u. The inset details the construction of Vtot 
with a short-range attractive potential Vsr and a long-range 
screened Coulomb repulsion Vy ■ 
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where Lb — e 2 / (AwEoekBT) is the Bjerrum length, £ = 
4:TtLb(pZo + 2p s )~ 1 / 2 is the Debye length, with p the pro- 
tein number density, Zq the net charge on a protein, p s 
the salt number density in the buffer and X a correction 
factor that depends on both <j> and £ (Eq. (11)-(15) in 
The resulting Yukawa potential gives a realistic de- 
scription of the effective repulsion between proteins for 
the relatively high volume fractions investigated. Here 
the explicit (^-dependence of X incorporates the effect 
of screening of a protein by other proteins, whereas the 
contribution of additional counter-ions with increasing p 
enters in the calculation of £. Under the chosen con- 
dition of low background electrolyte, we find it impor- 
tant to properly describe the change of £ and X with 
(j>. Under excess salt, or infinite dilution, the GOCM 
reduces to the repulsive part of the classical DLVO po- 
tential. Particles are assumed to interact simultaneously 
via Vtot = Vsr + Vy, shown in Fig.^for various <j> values. 
The integration time-step is fixed to At — 5 10 -3 in units 
of yjma 2 Je. A cutoff at r c = 8£ is applied to reduce the 
computational effort, without significantly altering the 
model, since V t ot(r c ) < 10~ 3 . 

Whilst the theoretical work of Belloni allows us to 
model the electrostatic repulsion without any fitting pa- 
rameter, no theory is capable to properly account for the 
attractive part. However, in the case of lysozyme and 
numerous other globular proteins, attraction is known 
to be short-ranged and of moderate strength [28l|. Pre- 
vious scattering studies on lysozyme suspensions under 
high salt conditions have shown that the attractive part 
of the effective potential does not significantly depend on 
temperature |13l l29l l30 | . Following these indications, the 
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FIG. 2: Comparison of the structure factors S(q) obtained by 
MD simulations (lines) with SAXS (symbols), for lysozyme 
suspensions of volume fraction 0.188 at corresponding tem- 
peratures for a — 90. 

depth and width of the Lennard-Jones potential in Eq. 
JTJ are chosen by running simulations at fixed cf> = 0.188 
for various a, and searching for the a-e values best match- 
ing the experimental peak positions of S(q) at 5°C. 

The results from simulations obtained with a = 90, 
corresponding to an attraction range of 3.6%<r and e = 
160meV or m 6.5fceT, are compared with the measured 
S(q) in Fig. Hat = 0.188. The position of the cluster 
peak and its T-dependence are properly accounted for the 
chosen a and e values. It is important to point out that 
also the location of the nearest neighbor peaks at « 2.25 
nm _1 coincides with previous SANS measurements Q. 
We also note that the low g-limit is rather well described, 
suggesting that the present potential properly accounts 
for the system compressibility (oc S(0) _1 ) for all temper- 
atures. This remarkable agreement infers that the use 
of an effective potential captures the essential features 
necessary to describe the clustering process in lysozyme 
suspensions. 

To validate the resulting potential, we turn to examine 
other densities. In Fig. |3 numerical results are compared 
with the corresponding experimental data at T = 5°C 
for volume fractions ranging from = 0.085 to 0.201. 
It is important to stress that the same well depth and 
range in the mixed potential is used at all investigated 
4>: the only (^-dependence in the potential arises from 
the repulsive part, where both £ and X vary (see FigQ. 
Again, the numerical S(q) properly reproduces the fea- 
tures observed experimentally. With increasing </>, the 
cluster peak position does not significantly change while 
its amplitude systematically decreases. Moreover, the 
agreement is quantitative for the cluster and nearest- 
neighbor peak positions as well as for their amplitudes, 
albeit slightly overestimated at the lowest (p. Also the 
compressibilities are well reproduced, with an initial de- 
crease with increasing </>, saturating to a roughly constant 
value. 

The ability of the potential to correctly reproduce the 
structure of the system suggests to analyze the equilib- 
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FIG. 3: (^-dependence of the static structure factor S(q) mea- 
sured by SAXS (upper panel) and obtained by MD simulation 
(lower panel) with potentials of FigQ The dependence of the 
Debye length £ on volume fraction is shown as inset. 

rium configurations provided by MD. Particles are con- 
sidered to be in the same cluster when the separation be- 
tween pairs of nearest-neighbor particles is smaller than 
r max — 1-065(7, a value corresponding to the distance 
at the maximum in V tot (see Fig^l. Fig. shows the 
resulting distribution n(s) of clusters of size s, for sev- 
eral state points. We note that these results are qual- 
itatively unaffected by different choices of the bonding 
distance. At the smallest <fi, the suspension is essen- 
tially composed of small clusters with a distribution that 
rapidly and monotonically decays to zero. As cj> increases, 
larger clusters form and the distribution progressively de- 
velops a power law behavior rc(s) oc s r , with an exponent 
t w —2.2 consistent with random percolation. 

We can try to exploit the micelle analogy in order to 
understand the strong variation of the cluster size distri- 
bution shown in FigQ] In the case of charged micelles 
at low ionic strength, one expects a dramatic change in 
the micellar growth scenario at a crossover concentration 
where the screening length becomes comparable to the 
micellar length [3lj . At this point, a sharp transition to 
accelerated micellar growth with a very broad size dis- 
tribution occurs. In our case, the crossover is reached 
once £/<7 w 0.5. For smaller £, we then observe enhanced 
cluster growth and a broadening of the size distribution. 

For (j> > 0.148 percolating states are found where the 
system forms a space-spanning cluster. From the analy- 
sis of the MD trajectories one observes that clusters are 
highly transient. Even above the percolation threshold, 
the system is still in a fluid phase, and percolation does 
not imply gelation, for which a long-living network is nec- 
essary. Interestingly, the observed transient percolating 
states resemble those of the transient networks found in 
the semidilute regime of charged micelles, where under 
the same salt-free conditions a persisting low-q peak is 
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FIG. 4: a) Cluster size distribution n(s) for various (f> at 
T — 5°C. The distribution for small sizes follows random 
percolation dependence of s~ 2 ' 2 . For <f> > 0.148 large clus- 
ters are found indicating that the system is percolating, b) 
Snapshot at <f> — 0.125 with few selected clusters of size 
n = 10,30,82,203. 

also present[32|. The observation of isolated clusters, as 
depicted in Fig^J)), reveals that the clusters are highly 
random. We find both space-filling spherical shapes as 
well as filamentous structures, made primarily of single 
chains, that are very different from those observed in 
other more repulsive systems [fil l23j . Moreover, in the in- 
vestigated region of phase space, there is not a preferred 
cluster size, since a peak in n(s) is never observed. 

Clusters observed under the present no-salt conditions 
are generated by the competition between attraction and 
repulsion, as clearly indicated by the presence of the pre- 
peak in the structure factor and should not be considered 
as the transient clusters that are commonly observed on 
approach to a phase-separation boundary in short-range 
potentials, whose spectral signature is a peak in S(q) 
at q — > 0. In the investigated T-range, attraction is 
always sufficiently counter-balanced by the electrostatic 
repulsion, allowing for the existence of stable clusters 
and preventing phase separation. The resulting clusters 
are highly polydisperse and their structural development 
with volume fraction leads to transient networking, with 
hallmarks of random percolation. 

In summary, we demonstrated that equilibrium cluster 
formation in complex protein solutions can be reproduced 
using a unifying " colloid-approach" , with a simple one- 
component effective potential between proteins composed 
of a short-range attraction and a longer range repulsion. 
This indicates the existence of a fundamental principle 
for self-assembly in biological solutions. 
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